Sequential Desynchronization in Networks of Spiking Neurons with Partial Reset 
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The response of a neuron to synaptic input strongly depends on whether or not it has just 
emitted a spike. We propose a neuron model that after spike emission exhibits a partial response to 
residual input charges and study its collective network dynamics analytically. We uncover a novel 
desynchronization mechanism that causes a sequential desynchronization transition: In globally 
coupled neurons an increase in the strength of the partial response induces a sequence of bifurcations 
from states with large clusters of synchronously firing neurons, through states with smaller clusters 
to completely asynchronous spiking. We briefly discuss key consequences of this mechanism for more 
general networks of biophysical neurons. 
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The brain processes information in networks of neu- 
rons, which interact by sending and receiving electrical 
pulses called action potentials or spikes. The response 
of a neuron to incoming signals strongly depends on 
whether or not it has just sent a spike itself. After the 
initiation of a spike the membrane potential at the cell 
body (soma) is reset towards some potential and the re- 
sponse to further synaptic input is reduced due to the 
refractoriness of the neuron The dendritic part of 
the neuron where incoming signals are integrated, is af- 
fected only indirectly by this reset due to intra-neuronal 
interactions [2, S U|- 

Several multi-compartment models have been pro- 
posed, in which different parts of a single neuron interact 
to characterize this effect [2|]. For instance, in a two- 
compartment model jlj of coupled dendrite and soma, 
the membrane potential at the soma is reset after spike 
emission while the dendritic dynamics is affected only 
by the resistive coupling from the soma to the dendrite. 
This accounts for the fact that in several kinds of neurons 
residual charge remains on the dendrite (following the so- 
matic reset), that is then transferred to the soma[3, 
Thus the dynamics of the individual neurons is modified 
which severely affects the collective capabilities of net- 
works of such neurons. 

In this Letter we propose a simple neuron model which 
captures the response to residual input charges following 
spike emission in form of a partial reset and at the same 
time allows an analytical study of the collective network 
dynamics. A fraction c G [0, 1] of the residual supra- 
threshold input charge is kept by the neuron after reset. 
For c = all additional input charge not needed to trigger 
a spike is lost after spike emission, whereas for c = 1 the 
total input charge is conserved 0. Although the regime 
< c < 1 is likely to be the biologically more relevant, 
to our knowledge, there are so far no systematic studies 
of the dynamics of networks of neurons with partial re- 
sponse. To reveal the basic mechanisms underlying the 
collective dynamics of networks of such neurons we fo- 



cus on networks of globally and homogeneously coupled 
neurons. Despite their simplicity these networks already 
exhibit a rich variety of dynamics that is controlled by the 
partial reset. In particular we find and show analytically 
that for a broad class of neurons there is a desynchro- 
nization transition in the network dynamics determined 
by a sequence of bifurcations: For small c the fully syn- 
chronous state coexists with a variety of cluster states 
(cf. 0), with differently sized groups of synchronously 
firing neurons. With increasing c, states with clusters of 
size a and larger become sequentially unstable at bifurca- 
tion points Cct satisfying < . . . < Cc? < < 1 such 

(2) 

that for sufficiently large c > C£ T , only an asynchronous 
state is left. We investigate the main mechanism generat- 
ing this sequence of bifurcations analytically and give an 
intuitive explanation. We also discuss key consequences 
of this novel desynchronization mechanisms for biophys- 
ically more detailed systems. 

Consider a network of N oscillatory neurons [3] , whose 
state at time t is characterized by a phase variable fa, 
i G {1,2, ... ,N}, that constantly increases with time 
4i(j)i = 1. The membrane potentials Ui — U ((/>i) are 
specified by a rise function U, that characterizes the 
subthreshold dynamics of a neuron. Here U is smooth, 
strictly monotonically increasing (U' > 0) and normal- 
ized to 17(0) = and £7(1) = 1. 

A neuron j generates a spike when its membrane po- 
tential crosses a threshold, Uj(t~) > 1, which in turn may 
trigger an avalanche of spikes (cf. Fig. [lj : Neurons reach- 
ing the threshold due to the free time evolution define the 
triggering set Oo = {j | Uj (t~) = 1}. The units j G Oo 
generate spikes which are instantaneously received by all 
the neurons i in the network. In response, their poten- 
tials are updated according to 
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where Eij > determines the strength of a directed 
synaptic connection from neuron j to i. The initial 
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pulse may trigger certain other neurons k G 0i = 
|fc | Ufe (t~) < 1 < to spike, etc. This process con- 

tinues n < N steps until no new neuron crosses the 
threshold. At each step m G {2,3,..., n} the potentials 
are updated according to 



(m+l) (m) 



(2) 



where m = w^"' ^ < 1 < u^ m H. The phases im- 
mediately after the avalanche 6 = Ug=o ®<? °f s i ze 



101 are obtained via 
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where R is the partial reset function. Here we focus on 
the linear form R (() = c£, with c G [0, 1] specifying 
the remaining fraction of supra-threshold input charges 
after reset. As a key example of the collective dynam- 
ics of neurons with partial reset, we here study neurons 
with convex rise function (U" > 0, modelling e.g. a class 
of conductance based integrate-and-fire neurons) , which 
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Figure 1: Model Dynamics, (a) Membrane potential u of a 
single neuron reaches threshold by intrinsic dynamics at time 
ti, sends a spike and is reset. In response to a sub-threshold 
input at t2, the potential is increased instantaneously by e. 
At time t% the neuron receives supra-threshold input, u (t^) + 
i > 1, sends a spike and a residual fraction c G [0, 1] of the 
charge, not used to reach the threshold, remains and yields 
a partial reset to the value u (t^ ) = c (u3 — 1). (b) Spike 
avalanche (a = 3, n = 3) in a network of N = 4 neurons, (i) 
The avalanche is triggered by neuron 1 reaching the threshold 
(defining the triggering set 0o = {1}) and sending spikes to 
all other neurons, (ii) This forces neuron 2 to cross threshold 
which then itself spikes (0i = {2}). (iii) In turn, this induces 
a firing of neuron 3 (03 = {3}), completing this avalanche, 
(iv) Finally the new states of all neurons are determined using 
the partial reset J3]) with O = {1, 2, 3}. 




Figure 2: Sequential desynchronization transition in a net- 
work of N = 50 neurons (17 = U b , b = -3, e = 0.0175). 
The phases (f>i of all neurons are plotted against the s-th 
spike of a reference neuron. Starting from a synchronous 
state and perturbing at s = 5, the phase dynamics are 
shown for (a) ci = 0.025 < c^r (inset: magnification), (b) 
c 2 = 0.5 G (c^'\cS ) ) and (c) c 3 = 0.7 > ciV . (d) Observed 

cluster sizes (dots) in the asymptotic dynamics of 6000 sim- 
ulations for each c value starting from different initial condi- 
tions, red line: exact theoretical prediction |TlO|| above which 
clusters are unstable. 



are homogeneously and globally coupled without self- 
interactions, Sij = (l-dij)e, and total input strength 
e = J2 j e ij = (N-l)i<l. 

Systematic numerical investigations indicate a strong 
dependence of the network dynamics on the partial 
reset strength c: In particular, we find synchronous 
states, cluster states, asynchronous states and a sequen- 
tial desynchronization of clusters with increasing c. More 
detailed, if c is sufficiently small, the long-term dynamics 
is dominated by many coexisting cluster states in which 
neurons are synchronized to differently sized groups re- 
sulting in a periodic state of the network (cf. Fig. [2]). As 
c increases, less and less clusters are observed with the 
maximal number of units per cluster decreasing. Above 
a critical value of c only an asynchronous state remains. 

What is the origin of this rich repertoire of dynamics 
and which mechanisms control the observed transition? 
To answer these questions, we analytically investigate the 
existence and stability of periodic states with clusters of 
arbitrary sizes. It turns out that the sequence of bifur- 
cations is controlled by two effects: sub-threshold inputs 
that are always synchronizing and supra-threshold inputs 
that may be synchronizing or desynchronizing depending 
on the strength c of the partial reset. 

As the first step we show that the fully asynchronous 
(splay 0) state exists and is stable for all c G [0, 1]. It 
is defined by identical inter-spike intervals between each 
pair of subsequently and individually firing neurons. A 
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firing map maps the phases 4>% of the system just before 
one avalanche to the phases just before the next. To con- 
struct this map for the asynchronous state we evaluate 
the current spike (a 1-neuron "avalanche") and shift all 
phases by the same amount a such that the largest of the 
resulting phases is at threshold. Without loss of general- 
ity, we label the neurons' phases in ascending order such 
that the phases <f>* and the shift a* uniquely define the 
asynchronous state; they are determined self-consistently 
by 4>\ = a* > and (\>\ = f/" 1 (U x ) + e) + a* for 
I e {2, . . . , N} such that (j)*N = 1- Homogeneity of the 
network implies invariance of such an asynchronous state 
under the firing map for every e < 1. Applying a small 



perturbation*^ ' = (s^, 
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to the N— 1 phases 



which are not at threshold and linearizing the firing map 



(cf. |lCllll]|) yields the perturbations after the next firing 
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Here A is a matrix whose only non-zero elements are 
A i+1 ,i = a, for i G {1, . . . N — 2} and A^-i = — ojv-i 
where 



u 1 m 



U> (U-i(U (#)+£)) 
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for % G {1, ... N - 1}. Since U' > 0, U" > 0, and e > 
we have Oj < 1. The Enestrom - Kakeya theorem [12J ap- 
plied to the matrix A implies that its eigenvalues satisfy 
| Af | < 1 for alH < iV— 1, showing that the asynchronous 
state is linearly stable. The stability properties of this 
state are identical for all c G [0, 1] because all neurons 
fire individually and do not initiate any avalanche with 
supra-threshold inputs. 

Next we investigate the stability properties of a peri- 
odic cluster state under the return map, i.e. the mapping 
of all phases just before the triggering of an avalanche 
to all phases just before the same avalanche reoccurs. 
Such a cluster state exists (i.e., is invariant) unless the 
maximal cluster size is too large such that this cluster ab- 
sorbs neurons not belonging to it or is absorbed by other 
clusters. Given that the specific neuron N G Oo stays 
in the triggering set of the avalanche, the return map M 
equals the map defined form the hyperplane 4>n = 1 to 
itself. It is fully specified by the number m of avalanches, 
1 < m < N, by the number a s , s G {1, . . . , to}, of 
neurons spiking in each avalanche, and by the subse- 
quent phase shifts a s that fix the time lags between the 
avalanches. This information is determined from the ini- 
tial phase vector (<t>i, . . . , 4>n-i, 1) and grouped into a 
firing sequence T = [(e a , u s )]^Lo' setting e s = a s e. For 
given T the return map then reads 



Mr (fa) = 5 CTm o H e 



■oS<j 2 oH £2 oS ai o J £l (fa) (6) 



for i £ 9. Here S a {<p) = <j> + & is the map mediat- 
ing a pure phase shift, H x ((f>) = U^ 1 (U ((f)) + x) speci- 
fies the sub-threshold response to an incoming spike and 



J x (4>) = U^ 1 (R (U (4>) + x — e — 1)) represents the par- 
tial response R to supra-threshold input. By definition 
we have Mj?(l) = 1. The conditions 

M r ([U- 1 (1 - ae),U~ l (1 - (a - l)g)]) 

C [U- 1 (1 - as) , 1] (7) 

for all a G {1, . . . , a{\ then ensure that all neurons firing 
in the first avalanche a\ will fire together in an avalanche 
after return of neuron N to threshold (4>n — 1)- Thus 
([7]) ensure stability of a cluster 6 of size |6| = a\. For 
general R and U these conditions yield upper and lower 
bounds [13] on the maximal size of a cluster to be stable 
under the return map. Here we focus on the specific rise 
function U b {4>) = \ In (l + (e b - l) cj)) , b < 0, where the 
change of phase differences due to sub-threshold inputs 
is independent of the phase, i.e. H e {4>) — H e (tp) = H £ o 
S a {4>) ~ H e o S a (ip) for a > 0. For i £ 9 the return map 
((6]) then simplifies to 



My ( 



with a = 1 — H, 



— So 



O H( N _ ai ) i o J aiS (fa) 



(8) 



> the conditions 
condition 



, Y_ai)g °Jaii (!)• Since M' F > and 
are all satisfied if the single 



M T (U^ 1 (1 - e)) > U^ 1 (1 - i) . 



(9) 



holds, where a single unit triggers the avalanche. A 
generic perturbation will disturb the cluster such that it 
gets triggered by a single unit only. Thus, if an avalanche 
of size a exist, condition ([9]) is sufficient and necessary for 
its stability. 

As a specific example, consider a linear partial reset 
function R(() — c(. Using equality in ([9]), the bifurca- 
tion values Cc?^ above which a cluster state with maximal 
cluster size a becomes unstable are then determined im- 
plicitly by the equation 



2 b(l-[(AT-a)+cW(a-l)]e) ^-be - l) = ( 



e-fc&'e _ i 



(10) 

Figure [2] shows an explicit example of these theoreti- 
cal predictions for the bifurcation values which well 
match the numerical results. 

In general, we infer from (flOf that 

o< c w< c (f- 1 )<...< c (?)<i (11) 

independent of specific parameters b, e and N . First, this 
implies that the entire sequence of bifurcations is guar- 
anteed to occur in the physically relevant open interval 
c G (0,1). Second, with increasing c, states with larger 
clusters become unstable before states with smaller clus- 
ters. In particular, the synchronous state becomes unsta- 
ble first and cluster states with at most two synchronized 
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Figure 3: Synchronizing and desynchronizing mechanisms 
underlying the desynchronization transition. The phase- 
potential relation are shown as solid black curves. The phase 
differences before (A - ) and after spike reception (A + ) are 
indicated, (a) Phase synchronization due to sub-threshold 
inputs coactig with either (b) synchronization for strongly 
refractory partial resets (c w 0) or (c) desynchronization for 
conservative partial resets (c « 1) determine the stability of 
clusters under the return map (J6j> - 



neurons become unstable last. Third, for a = 2 we find 
that the largest bifurcation point 



.(2) _ 



X 1 



o b(l-(N-l)e) 



(12) 



can be arbitrarily small, e.g. as b — * — oo. Thus the 
entire sequence of desynchronizing bifurcations can occur 
for arbitrary small c. 

The from Q of the return map reveals the mecha- 
nisms underlying the desynchronization transition as the 
interplay between synchronization due to sub-threshold 
inputs mediated by -ff(jv- a )£ (cf. Fig. ((3j)a) and fur- 
ther synchronization or desynchronization due to supra- 
threshold inputs and partial reset mediated by J a g, de- 
pending on the strength of the partial reset (cf. Fig. 
([3])b,c). The large clusters get unstable first since they 
receive less synchronizing sub-threshold inputs from the 
other smaller clusters and additionally the desynchro- 
nization due to the reset is stronger in larger avalanches. 

The observed desynchronization transition prevails for 
networks of inhomogeneously coupled units and neurons 
with rise functions of mixed convex and concave curva- 
ture, as e.g. characteristic for quadratic integrate-and- 
fire neurons the normal form of type I excitable 

neurons. Moreover, our simple model system can be con- 
nected to biophysically more detailed type I models by 
comparing spike time response curves that encode the 
shortening of the inter-spike intervals (ISI) following an 
excitatory input at different phases of the neural oscilla- 
tion. An excitatory stimulus that causes the neuron to 
spike will maximally shorten the ISI in which the stimu- 
lus is applied. Additionally the following ISI is typically 
affected as well. This effect can be characterized by an 
appropriately chosen partial reset in our simple system 
[101 ] - Networks of two-compartment conductance based 
neurons indeed exhibit similar desynchronization transi- 
tions when varying the coupling between soma and den- 
drite (not shown) which in our simplified model controls 



the partial reset. 

In summary, we introduced a simple model of spik- 
ing neurons with partial reset to investigate collective 
network effects of possible residual charges that may be 
important after somatic reset. Already for globally and 
homogeneously coupled networks we find that residual 
charges present after spike generation drastically affect 
the network dynamics. We revealed a new desynchro- 
nization mechanism that controls a sequential destabi- 
lization of cluster states. For no or only small fractions 
c 6 0, Cc^^j of conserved charge, the synchronous state 
and cluster states with many different cluster sizes coexist 

Ccr , 1 , only the asyn- 
chronous state is left. For intermediate c £ c^r jCcr 
there is a sequence of bifurcations, each destabilizing the 
largest stable cluster. Interestingly, this entire sequence 
may occur in an interval at arbitrarily small c > 0. 

The mechanism for neural desynchronization discussed 
above differs strongly from known mechanisms that are 
based, e.g., on heterogeneity, noise, or delayed feedback 
[lR fig . Possibly, this novel mechanism may also be used 
in modified form to prevent synchronization in neural 
systems like in Parkinson tremor or in epileptic seizures 
[lf|]. This calls for a future systematic study of the im- 
pact of c and related parameters that modulate local re- 
sponse properties and thus synchronization. The simple 
model system introduced above offers the first example 
of an analytically tractable network model which, based 
on partial reset, characterizes an essential feature of bio- 
physically detailed compartmental models 0, [TtJ • 
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